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Abstract 


The  implications  of  Segletes’  recent  equation  of  state  are  examined,  when  idealized  m  the 
condition  defined  as  quasi-harmonicity.  Results  indicate  that  at  large  volumetric  strains,  a 
proportionality  no  longer  holds,  in  general,  between  the  volumetric  and  vibrational  stiffnesses. 
The  governing  relation  between  these  two  stiffnesses  is  presented  and  is  a  function  of  the 
characteristic  frequency  of  the  lattice,  alternately  expressed  in  terms  of  the  lattice  spacing.  It  is 
further  shown  that  the  quasi-harmonic  idealization  of  Segletes’  equation  will  approach  the 
harmonic  approximation  in  the  limit. 
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1.  Background 


The  harmonic  approximation  has  been  the  traditional  method  to  analyze  atomic  vibrations  and 
lattice  behavior,  dating  back  to  Planck.  Several  years  later,  in  the  early  1900s,  Einstein,  Debye, 
and  Bom  and  von  Karman  [1]  used  the  approximation  to  develop  theories  of  lattice  specific  heat. 
In  the  years  that  followed,  Brillouin  [2]  also  employed  the  harmonic  approximation  to  analyze 
wave  propagation  in  lattice  structures. 

In  the  harmonic  approximation,  atoms  vibrate  with  sinusoidal  motion.  In  order  for  an  atom 
within  the  crystal  lattice  to  retain  this  harmonic  behavior,  the  specific  potential  energy  well  of 
the  lattice,  Ep^lr),  must  be  parabolic  in  shape,  where  r  is  the  coordinate  of  the  atom  in  question 
within  the  lattice.  SpecificaUy,  EEp„,~(r-r^f.  Under  these  conditions,  the  restoring  force  to 
the  atom  is  proportional  to  its  displacement  from  the  rest  position,  and  the  resulting  oscillatory 
motion  of  the  atom  is  sinusoidal — ^that  is  to  say,  harmonic.  The  harmonic  motion  of  an  atom  has 
a  single  frequency  associated  with  it,  iOq,  and  that  frequency  is  calculable,  just  as  in  macroscopic 

vibration  theory,  by  the  relation  (Og  =  (-3P/3r)  /  M ,  where  M  is  the  mass  of  the  atom,  P  is  the 

restoring  force,  and  -BP/Br  is  the  local  spring  constant  in  the  direction  of  the  disturbance,  which, 
for  the  harmonic  case,  is  independent  of  position.  If  a  lattice  of  vibrating  atoms  is  considered, 
rather  than  just  a  single  oscillating  atom,  a  whole  spectrum  of  frequencies  is  shown  to  arise  in 
the  lattice.  Furthermore,  in  a  three-dimensional  (3-D)  lattice,  for  every  frequency  of  the  spectrum 
associated  with  a  disturbance  in  the  direction  of  the  wave  propagation  (a  longitudinal  wave),  there 
are  two  vibrational  components  associated  with  disturbances  perpendicular  to  the  direction  of 
wave  propagation  (shear  waves). 

Let  us  remember,  now  and  throughout  this  report,  that  vibrations  are  the  means  by  which 
internal  (i.e.,  potential)  and  kinetic  energy  are  exchanged  in  a  spring.  Every  type  of  vibration, 
no  matter  how  complicated  kinematically,  has  an  associated  spring  constant  (or  alternately,  wave 
speed)  and  displacement  variable.  The  derivative  of  the  spring’s  potential  energy  with  respect 
to  the  displacement  variable  is  known  as  the  “force,”  and  the  derivative  of  the  force  with  respect 
to  this  displacement  variable  is  the  local  spring  constant  (or  stiffness),  which  governs  the 
vibrational  frequency  for  infinitesimal  vibrations  of  the  spring. 
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An  early,  important  result  of  lattice-vibration  theory  was  in  showing  that  vibrational  energy 
in  a  lattice  can  manifest  itself  as  thermal  pressure  only  if  the  lattice  vibrational  frequencies 
change  with  lattice  deformation — ^that  is  to  say,  if  the  potential  energy  well  is  not  parabolic  in 
shape.  This  condition  of  nonparabolicity  is  referred  to  as  anharmonicity.  Likewise,  nonparabolic 
terms  in  a  Taylor  expansion  of  the  potential  energy  well  about  the  reference  lattice  spacing  are 
called  anharmonic  terms.  In  practical  terms,  the  empirical  metric  for  anharmonicity  usually 
entails  measuring  the  nonparabolicity  of  the  lattice’s  energy  well  only  at  the  reference  lattice 
spacing  (e.g.,  the  method  used  by  Rose  et  al.  [3]).  Here  and  throughout  this  report,  the  reference 
condition  is  one  of  zero  temperature  and  pressure. 


Quantitatively,  a  lattice’s  vibrational  energy  gets  coupled  into  thermal  pressure  by  way  of  the 
Griineisen  function,  macroscopically  defined  as  F  =  V(dpldE)y,  where  V is  the  material’s  specific 

volume,  and  p  and  E  are  the  pressure  and  specific  energy  of  the  material,  respectively.  From 
lattice  theory,  there  is  a  Griineisen  value  associated  with  each  vibrational  frequency  in  the  lattice, 
given  by  the  following  expression  and  taken  at  constant  temperature,  T: 


fa.,1 


JT 


(1) 


where  F,  is  the  Griineisen  value  associated  with  the  frequency,  (0„  of  a  longitudinal  or  transverse 
(i.e.,  shear)  vibration.  Were  the  F,-  values  for  each  frequency  truly  independent  of  each  other,  the 
macroscopically  observed  F  value,  being  an  aggregate  of  all  the  F,-,  would  be  very  difficult  to 
characterize  in  terms  of  the  individual  F,.  However,  one  fortunate  simplification  to  the  problem 
arises  from  Griineisen’s  assumption,  which  was  originally  made  based  on  empirical  data  and  later 
supported  theoretically  by  the  lattice  theory  of  Debye.  The  assumption  asserts  that  the  whole 
spectrum  of  lattice  vibrational  frequencies  changes  in  a  characteristically  similar  fashion  with 
lattice  spacing,  so  that  the  knowledge  of  how  one  characteristic  frequency  changes  with  volume 
defines  the  behavior  of  the  complete  vibrational  spectrum.  The  effect  on  eqn  (1)  is  to  remove 
the  ‘i’  subscript,  with  the  resulting  equation  being  characteristic  of  the  behavior  of  the  complete 
vibrational  spectrum.  From  the  early  theories  of  specific  heats,  this  characteristic  frequency  (a 
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microscopic  quantity)  was  shown  to  be  directly  proportional  to  a  characteristic  temperature,  0 
(a  macroscopic  quantity).  Thus,  it  is  common  to  find,  in  the  literature,  an  interchangeable  use 
of  the  expressions  ((o/cOq)  ^nd  (©/0o),  where  the  “0”  subscript  refers  to  the  reference  value. 
Another  simplification  is  afforded  by  the  approximation,  supported  by  data,  that,  for  solids  in  the 
absence  of  phase  change,  eqn  (1)  is  independent  of  temperature,  so  that  the  partial  derivatives 
in  eqn  (1)  become  ordinary.  In  this  case,  T,  (o,  and  0  all  become  functions  of  specific  volume 
only.  Work  in  the  area  of  a  temperature-dependent  Griineisen  function  was  done  by  both 
Blackman  [4]  and  Barron  [5],  who  examined  and  modeled  the  variation  of  F  with  temperature 
at  low  values  of  temperature  (r/0  «  1)  and  found  mild  variations.  Later,  Grodzka  [6] 
summarized  work,  using  the  compression  of  porous  media,  which  indicated  that  some  materials 
(notably  aluminum)  seem  to  have  a  significant  decrease  in  F  at  elevated  temperatures.  Aluminum 
aside,  the  temperature  independence  of  the  Griineisen  parameter,  while  empirical,  seems  valid 
over  a  wide  range  of  conditions  and  solids. 

Even  with  the  simplifications  afforded  by  the  temperature-independent,  coupled  frequency 
form  of  eqn  (1),  the  difficulty  still  arises  of  how  to  describe  the  anharmonic  behavior  of  lattices, 
given  that  the  prevailing  theories  describing  their  fundamental  behavior  are  based  on  harmonic 
approximations.  Slater  [7]  and  later  Dugdale  and  MacDonald  [8]  used  the  harmonic 
approximation  as  a  basis  to  develop  anharmonic  theories  of  the  Griineisen  function. 

To  better  understand  the  model  of  Dugdale  and  MacDonald  [8],  as  well  as  more  recent  work, 
let  us  construct  a  framework  for  analyzing  anharmonic  vibrations.  Given  a  simple  spring-mass 
system,  consider  replacing  the  linear  (harmonic)  reference  spring  with  a  nonlinear  (i.e., 
anharmonic)  spring.  Vibration  theory  may  be  used  to  compare  the  frequencies  of  vibration  for 
these  two  systems.  Though  the  frequency  of  vibration  for  a  nonlinear  spring  changes  with 
amplitude,  the  situation  of  low-temperature  atomic  vibrations  may  be  likened  to  the  situation  of 
a  spring  experiencing  infinitesimal  vibrational  amplitudes  superimposed  about  a  nominal  spring 
extension.  In  fact,  the  notion  of  thermal  expansion  in  a  solid  is  directly  analogous  to  a  nonlinear 
spring,  in  that,  as  the  temperature  is  raised,  the  vibrational  amplitude  increases,  and  as  the 
amplitude  increases,  the  nonlinearity  of  the  spring  causes  the  average  spring  extension 
(i.e.,  density)  to  change,  even  though  the  nominal  applied  force  remains  constant.  The  nominal 
spring  extension  is  analogous,  in  a  lattice,  to  the  atomic  spacing  relative  to  the  reference 
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configuration,  and  the  small  vibrational  oscillations  that  occur  do  so  about  this  deformed 
configuration.  Under  this  condition  of  infinitesimal  vibrations,  corresponding  to  the  low- 
temperature  condition  of  the  lattice,  the  vibrational  frequency  of  a  nonlinear  spring  is  obtainable 
by  way  of  perturbation  theory  (i.e.,  by  linearizing  the  problem  and  considering  small 
perturbations  in  displacement  about  the  nominal  extension).  For  a  given  mass,  the  local  spring 
stiffness  determines  the  vibrational  frequency  of  the  spring  and  is  given  by  the  local  slope  of  the 
force-displacement  curve  for  the  spring,  which  is  a  function  only  of  the  current  spring  extension. 
Thus,  the  relationship  that  governs  how  the  frequency  response  is  altered  by  substituting  the 
anharmonic  spring  for  the  harmonic  reference  spring  is  given  as 


1 

M 


(2) 


where  the  “re/”  subscript  refers  to  the  reference,  harmonic  spring. 

To  further  extend  the  spring  analogy  to  a  lattice  of  atoms,  consider  the  case  of  a  3-D  lattice 
of  spring-connected  identical  masses.  In  the  general  case,  the  lattice  may  support  physical 
motion  both  in  the  direction  of,  as  well  as  transverse  to,  the  direction  of  wave  propagation. 
Furthermore,  to  be  general,  we  allow  the  lattice/spring  topology  to  include  diagonally  connected 
masses,  such  that  a  displacement  in  one  direction  may  induce  force  and/or  motion  in  the 
perpendicular  directions.  In  this  general  case,  the  vibrational  frequency  spectrum  of  the  spring- 
connected  lattice  includes  components  due  to  both  longitudinal  and  transverse  vibrations.  Such 
a  frequency  spectrum,  aggregated  of  both  the  longitudinal  and  transverse  spectra,  is  directly 
analogous  to  that  found  in  actual  solids.  Brillouin  [2]  has  shown  that,  if  a  Debye  frequency 
spectrum  is  assumed  for  a  material,  the  aggregated  characteristic  temperature  (frequency)  is 
describable  as  a  composite  of  the  longitudinal  and  transverse  frequency  spectra  as  follows. 

J_  +  A  =  J_  .  (3) 

0;  0^ 
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Others  [9, 10]  have  accomplished  a  similar  result  by  relating  an  aggregate  of  longitudinal  and 
transverse  wave-propagation  velocities  (related  to  their  associated  frequency,  to,  and  wavelength, 
X,  via  C  =  toX)  directly  to  the  Griineisen  function  (related  to  the  frequency  spectrum  by  eqn  [1]). 
For  example,  Pastine  [9]  provided  this  relation  as 

r  =  1/3  [1/3  -  (V/q)  dC/dV]  +  2l3[ll3-(yiC^)  dCJdV]  ,  (4) 

where  Q  and  C,  refer  to  the  longitudinal  and  transverse  wave  speeds,  respectively.  By  relating 
these  component  wave  speeds  to  variations  in  the  Poisson  ratio,  Pastine  was  able  to  derive  a 
Griineisen  expression  as  a  function  of  Poisson’s  ratio  and  showed,  for  the  special  case  of  a 
monatomic  face-centered  cubic  crystal,  a  result  quite  close  to  the  theoretically  derived  result  of 
Barron  [5].  However,  for  the  general  case,  no  guidance  was  offered  on  appropriate  selection  of 
functional  forms.  Vashchenko  and  Zubarev  [10]  adopted  a  similar  approach,  but  assumed  a 
particular  fitting  form  for  the  Poisson  ratio,  in  terms  of  the  cold  pressure  curve  and  a  fitting 
parameter. 

Unlike  these  other  approaches,  Plendl  [11]  considered  the  measured  frequency  spectra  for  real 
materials,  as  opposed  to  using  idealized  relationships  like  eqns  (3)  and  (4).  Rather  than 
attempting  a  decomposition  of  the  spectrum  into  longitudinal  and  transverse  components,  he  was 
able  to  bypass  such  an  approach  and  show  that  the  characteristic  frequency  of  a  material  is  that 
frequency  at  the  center  of  gravity  of  the  frequency  spectrum,  /((o),  the  so-called  “centro- 
frequency,”  which  is  given  by 


W2 

jo)  /((o)cfco 
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where  the  frequency  range  of  integration  is  selected  to  comprise  the  entire  range  of  the  acoustical 
vibration  spectrum. 

The  current  analysis  follows  the  lead  of  Plendl,  in  that  no  explicit  decomposition  of  the 
vibrational  spectrum  into  longitudinal  and  transverse  components  will  be  performed.  The  mere 
fact  that  a  characteristic  frequency  may  be  aggregated,  in  some  manner,  from  the  longitudinal  and 
transverse  frequency  spectra  (by  way  of  eqn  [3],  [4],  [5],  or  otherwise)  allows  us  to  introduce  the 
concept  of  a  vibrational  spring  constant,  which  is,  by  definition,  the  spring  constant  that 

characterizes  the  aggregated  (longitudinal  plus  transverse)  vibrational  frequency  spectrum.  If  a 
lattice  had  no  diagonal  coupling  and  didn’t  support  transverse  (shear)  waves,  the  characteristic 
vibrational  stiffness  would  simply  reduce  to  the  stiffness  of  the  component  longitudinal  spring. 
The  vibrational  spring  constant  will,  like  the  stifftiess  of  a  nonlinear  spring,  be  a  function  of  the 
nominal  lattice  spacing.  By  way  of  this  construct  of  the  vibrational  spring  constant,  eqn  (2)  may 
be  generalized  to  case  of  a  3-D  coupled  lattice: 


tO^  _  ^re/ 


(6) 


where,  now,  o)  is  a  characteristic  frequency  of  the  vibrational  spectrum,  rather  than  the  specific 
frequency  of  the  component  spring. 

Let  us  analyze  a  lattice  where  the  specific  unit  lattice  spacings  (i.e.,  length  per  cube-root  unit 
mass)  along  the  principal  lattice  directions,  x,  y,  and  z,  are  expressed  in  terms  of  the  specific 
volume  of  the  unit  cell  as  V=‘k^  =  xyz,  where  X  is  an  averaged  unit  lattice  spacing  per  cube-root 
unit  mass.  For  the  case  of  zero-temperature  volumetric  compression  or  distention,  a  specific 
force  may  be  defined  as  that  force  (per  2/3  power  unit  mass)  acting  on  the  face  of  a  unit  cell  of 
the  lattice,  given  as  definition,  p,,  is  the  cold  (0°  isotherm)  pressure  and, 

thus,  the  force  is  positive  in  compression. 

If  nothing  else  were  known  of  the  vibrational  frequency  spectrum,  it  would  be  natural  to 
investigate  the  proportionality  of  dF^,,ldx  and  dFIdX,  since,  for  a  one-dimensional  (1-D)  harmonic 
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lattice,  these  two  quantities  would,  in  fact,  be  identically  equal.  The  assumption  that  the 
characteristic  vibrational  stiffness  of  a  3-D  spring  can  be  made  proportional  to  an  associated 
volumetric  stiffness  has  been  pervasive  throughout  the  literature.  When  Debye  theory  indicated 
that  the  characteristic  frequency  varied  as  an  aggregate  of  the  elastic  wave  speeds.  Slater  [7] 
assumed  a  constant  Poisson  ratio  and  related  the  aggregated  vibrational  stiffness  (in  terms  of 
wave  speeds)  directly  to  volumetric  compressibility  (although  Slater  actually  concluded  something 
closer  to  [but  not  quite]  — dpIdX  rather  than  -dFIdX).  Brillouin  [2]  derives  results  utilizing 

(ii^-d^EJdX^,  the  right-hand  side  being  merely  -dF/dk.  It  will  be  shown  that  Dugdale  and 

MacDonald  [8]  also  utilized  this  proportionality.  More  recently,  Guinea  et  al.  [12]  explicitly 
assumed  this  proportionality  when  trying  to  infer  a  relationship  between  the  universal  cold  curve 
and  the  Griineisen  function. 

The  effect  of  such  an  assumption  of  proportionality  on  eqn  (6)  is  that  a  dF/dk  term  may  be 
substituted  for  each  dF/dx.  With  this  substitution,  the  derivatives  with  respect  to  k  are 
convertible  to  volume  derivatives  through  the  chain  mle, 

dF/dk  =  dF/dV -dV/dk  =  ,  (7) 

where  the  prime  symbol  (' )  denotes  ordinary  differentiation  with  respect  to  specific  volume. 
With  this  substitution,  eqn  (6)  becomes 

2 

CO^  _  ('g') 

Tpy^ "  Tp^  ■ 


If  one  takes  the  volume  derivative  of  eqn  (8)  and  divides  the  result  by  eqn  (8)  itself,  one  obtains 


^  dm/dV 
(0 


(py^^y '  ' 


(9) 
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The  first  term  from  each  side  of  the  equation  is  related  to  the  Griineisen  function  and  may  be 
substituted,  using  the  temperature-independent,  characteristic-frequency  form  of  eqn  (1),  to  obtain 


^r_  iprr  , 

V  V  ■ 


At  this  point,  the  values  for  the  harmonic,  reference  spring  may  be  substituted.  As  already 
mentioned,  the  value  for  F  in  a  harmonic  system  is  identically  zero,  since  O)  is  constant,  and  the 
last  term  involving  the  pressure  behavior  of  the  reference  system  may  be  acquired  using  the 
definition  p^^  =  -dE^^ldV,  where  is  the  parabolic  energy  well  expressed  in  terms  of  V  rather 
than  X.  The  result  is  the  well-known  Dugdale-MacDonald  relation: 


r 


1  _  viPcV^r 

3  2  (py^y  ■ 


(11) 


Inserting  the  harmonic  energy  potential  into  this  equation  results  in  a  Griineisen  value  of  zero. 
Such  a  result  illustrates  the  criticism,  by  Dugdale  and  MacDonald  [8],  of  the  model  of 
Slater  [7] — ^namely,  that  inserting  the  harmonic  cold  curve  into  Slater’s  model  does  not  produce 
the  required  value  of  F  =  0. 

This  report  will  show,  on  the  basis  of  Segletes’  recently  proposed  equation  of  state,  that,  for 
a  3-D  lattice,  the  presumed  proportionality  between  the  volumetric  and  vibrational  stiffness, 
which  has  pervaded  the  literature  and  leads  to  the  model  of  Dugdale  and  MacDonald,  is  strictly 
true  only  in  the  harmonic  limit  of  approximation.  Force-based,  vibrational  and  volumetric  moduli 
will  be  compared,  for  the  idealized  case  of  quasi-harmonicity  (defined  in  the  next  section),  with 
the  more  general  case  to  be  addressed  in  a  subsequent  report.  Finally,  the  quasi-harmonic 
idealization  of  Segletes’  equation  of  state  will  be  shown  to  approach  the  harmonic  approximation 
as  the  Griineisen  function  approaches  zero  in  the  limit. 
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2.  Segletes’  Equation  of  State  and  the  Quasi-Harmonic 
Idealization 

Recently,  Segletes  [13,  14]  proposed  an  equation  of  state,  of  Griineisen  form,  that  properly 
captures  both  the  lattice-potential  and  thermal  pressure  behavior  for  a  variety  of  crystalline 
materials,  into  the  megabar  range  of  pressures.  The  model  was  inspired  by  the  universal  cold 
curve  of  Rose  et  al.  [3].  However,  unlike  the  model  of  Rose  and  other  recent  equation-of-state 
models,  {e.g.,  Vinet  et  al.  [15,  16];  Baonza,  Caceres,  and  Ndnez  [17];  and  Baonza  et  al.  [18]), 
the  model  of  Segletes  casts  the  lattice  energy  potential  in  terms  of  the  characteristic  temperature 
of  the  lattice  {i.e.,  the  lattice  vibrational  frequency),  rather  than  lattice  spacing  (density).  The 
equation  of  state  is  given  as 


py\f-E  =  £:J[(0/©o)*^-  l]  +  ln(©/0o)} 


(12) 


In  this  equation,  p  and  E  are  the  pressure  and  specific  internal  energy  respectively,  K  is  a 

parameter,  given  by  where  Q  is  the  reference  bulk  sound  speed  at  zero 

temperature  and  pressure,  Ej,  is  the  specific  lattice  binding  energy,  and  tj/'  =  WF  =  (dEldp)y  is  a 
thermodynamic  variable  introduced  by  Segletes  [19,  20]  for  ease  in  manipulating  the  governing 
equations.  In  terms  of  this  variable,  tj/',  the  temperature-independent,  coupled-frequency  form  of 
eqn  (1)  is  given  as 


0  \|r 


(13) 


where  the  proportionality  between  ©  and  to  has  been  utilized.  The  generalized  cold  curve 
associated  with  eqn  (12)  is 

E^  =  E^a  -  [1  -  /s:in(©/©o)]}  (14) 
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and 


n  =  :^(0/0„f  ln(0/0o)  .  (15) 

\|/ 

When  the  K  parameter  takes  on  a  value  identically  equal  to  unity  and  the  Griineisen  function 
varies  as  F  -  (i.e.,  F-X),  the  equation  takes  on  behavior  that  will  hereafter  be  called  quasi- 

harmonic.  If  we  denote  the  variables  of  the  quasi-harmonic  equation  of  state  as  hatted 
quantities  the  quasi-harmonic  cold  curve  becomes 

[1  -  ln(^/<§»o)]}  (16) 


with  the  Griineisen  requirement  translating  to 


~  .  (17) 

Eqns  (16)  and  (17)  are  interpreted  as  quasi-harmonic  for  a  niunber  of  reasons.  It  is 
immediately  obvious  that,  when  K  equals  unity,  many  terms  drop  out  from  the  general  equation 
of  state,  eqn  (12),  and  the  equations  become  greatly  simplified.  Furthermore,  Segletes  and 
Walters  [21,  22]  examined  the  equation  of  state  using  a  simple  power  law,  V',  to  model  the  \|r 
function.  They  noted  that,  over  a  wide  range  of  anharmonicity,  the  exponent,  x,  required  to  fit 
well  the  model  of  Rose  et  al  [3],  varied  linearly  with  Rose’s  anharmonicity  parameter,  T). 
Though  it  did  not  match  precisely  in  the  limit,  an  exponential  relation  of  t);  ~  V  was  the  limiting 
trend  for  the  case  where  Rose’s  anharmonicity,  t|,  approached  zero.  Segletes  and  Walters  also 

showed,  for  the  case  of  1  and  tjlr  ~  7^^  (the  quasi-harmonic  case),  that  the  Segletes  equation 
of  state  satisfies  the  relation 

(py2/3)  +  [2  +  2/3(\ff/7)]i$r(py2/3)'  +  ^(^py7J^yl  =  q  ,  (18) 

whereas  the  Dugdale-MacDonald  relation,  eqn  (11),  noted  for  its  harmonic  limiting  behavior. 
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satisfies 


[2 + 113  (tif/v )]  v  {pv^^y  +  {pv^y  =  0  .  (19) 

Eqns  (18)  and  (19)  differ  by  only  the  first  term  of  eqn  (18). 

In  a  subsequent  report,  Segletes  will  treat  the  more  general  form  of  his  equation  of  state,  by 
showing  that  the  K  parameter  in  eqns  (12),  (14),  and  (15)  must,  for  real  materials,  be  a  (slowly 
varying)  function  of  volume,  which  asymptotes  to  the  quasi-harmonic  case  in  the  large-volume 
limit.  In  the  meantime,  the  cause  of  the  disparity  between  the  idealized  form  of  Segletes’  model 
and  that  of  Dugdale  and  MacDonald  [8]  will  be  shown  to  arise  firom  the  fact  that,  for  large 
volumetric  strains,  the  volumetric  and  vibrational  stiffiiesses,  dF/dX  and  /dx,  respectively, 
are  not  proportional. 

3.  Force-Based  Moduli 

The  interactions  within  a  lattice  are  governed  by  forces  acting  over  distances.  On  this  basis 
alone,  force  and  distance  would  seem  to  be  the  natural  variables  governing  the  interactions  within 
a  lattice.  More  compelling,  however,  is  the  fact  that  Segletes’  equation  of  state  shows  a  very 
distinct  preference  for  being  analyzed  in  a  force-based  framework.  Eqn  (18)  is  just  one  example 
of  how,  timft  and  again,  the  term  grouping  (PcV^\  which  is  the  force  on  the  face  of  the  lattice 
unit  ceU  (per  2/3  power  unit  mass),  appears.  It  is  for  this  reason  that  subsequent  results  will  be 
analyzed  in  a  force-based  framework,  by  which  it  is  meant  that  moduli  will  be  expressed  in  terms 
of  force  gradients,  rather  than  stress  gradients.  To  this  end,  force-based  moduli  will  be  defined 
for  a  zero-temperature  lattice,  at  an  arbitrary  volumetric  compression  or  distention,  under  the 
influence  of  an  infinitesimal,  superimposed,  elastic  disturbance. 

In  a  force-based  frame,  a  modulus,  with  units  of  stress,  may  be  associated  with  a  given  spring 
constant.  From  considerations  of  dimensional  analysis,  the  modulus  should  be  of  the  form  spring 
constant  divided  by  characteristic  length.  Alternately,  a  force-based  modulus,  associated  with  the 
volumetric  spring  constant,  dFIdX,  for  example,  may  be  defined  by  starting  with  the  stress-based 
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modulus  definition  (i.e.,  bulk  modulus),  substituting  a  force  term  for  the  stress  (pressure)  term, 
and  dividing  the  result  by  the  specific  area  (i.e.,  area  per  2/3  power  mass)  over  which  the  force 
is  acting.  For  a  lattice  that  is  only  deformed  volumetrically,  this  area  is  the  specific  volume 
to  the  2/3  power.  Thus, 


1  dF  ^  _  1 

dVIV  dV  3X  dX  ' 


(20) 


Comparison  of  the  stress-based  bulk  modulus,  B,  and  the  force-based  one  in  eqn  (20)  reveals  that, 
in  general,  Bp=B-  2pJ3.  Under  conditions  where  the  pressure  and  temperature  are  zero,  the 
force-based  volumetric  modulus  will  equal  the  stress-based  bulk  modulus  in  value. 

Turning  now  to  the  vibrational  spring  constant,  dF^ldx,  an  associated  force-based  modulus, 
call  it  Jp,  may  be  defined.  As  in  the  case  of  the  vibrational  spring  constant  itself,  this  modulus 
implicitly  includes  the  effects  of  both  longitudinal  and  shear  waves.  For  a  cubic  lattice,  it  is 
given,  in  terms  of  X,  by 


f  1 


\^yz  dx/x 


Jx^y^z^X 


1 ^ 
X  dx 


(21) 


4.  The  Relation  Between  the  Volumetric  and  Vibrational  Moduli 

We  will  summarize  what  is  known  of  the  problem  at  hand  by  enumerating  a  list  of  statements 
already  discussed,  so  that  subsequent  reference  may  conveniently  be  made  to  the  ideas  contained 
in  the  following  statements: 

(1)  In  a  Griineisen  material: 

(a)  the  vibrational  spectrum  of  a  lattice  is  composed  of  longitudinal  and  transverse  (shear) 
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vibrational  components.  This  spectrum  may  be  characterized  either  in  terms  of 
characteristic  wave  speeds  or  by  a  characteristic  frequency; 

(b)  the  characteristic  frequency  of  the  vibrational  spectrum  varies  with  volume  in  the  same 
manner  as  the  characteristic  temperature,  (O  ~  0;  and 

(c)  the  characteristic  temperature,  0,  is  related  to  the  \|;  function  (and  thus  F)  by  way  of 
eqn  (13). 

(2)  The  characteristic  frequency,  (O,  of  the  vibrational  spectrum  of  a  lattice  may  be  directly 
related  to  an  associated  spring  constant,  designated  dF^i,/dx,  of  a  material,  such  that  the 
characteristic  frequency  varies  as  ~  /dx.  This  spring  constant  may  be  alternately 
expressed  as  a  vibrational  modulus,  designated  Jf,  and  defined  in  eqn  (21). 

(3)  The  volumetric  stiffness,  dF/dk,  is  directly  derivable  from  a  material’s  equation  of  state,  and 
may  be  alternately  expressed  in  terms  of  a  volumetric  modulus.  Bp,  defined  in  eqn  (20). 

(4)  Segletes’  equation  of  state,  which  fits  the  behavior  of  crystalline  solids  very  well,  expresses 
the  behavior  of  the  lattice  in  terms  of  relative  characteristic  temperature,  0/0o.  Since  the 
characteristic  temperature  (frequency)  function  is  governed  by  the  vibrational  stiffness, 
dF^jjjIdx,  this  equation  of  state  may  be  used  to  compare  directly  the  volumetric  and  vibrational 
stiffnesses. 

(5)  Segletes’  model,  eqn  (12),  may  be  idealized  to  an  important  special  case,  termed  the  quasi¬ 
harmonic  idealization,  which  is  characterized  by: 

(a)  a  value  of  the  parameter,  K,  equal  to  unity;  and 

(b)  a  function  that  varies  as 

From  the  cold-pressure  curve  of  Segletes’  equation  of  state,  eqn  (15),  derivatives  may  be 
taken,  to  show  that 
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E 

=  --^—2 — 
\jr 


(@/0of  {[K+\\y-2/3(y\f/V)]  ln(©/0o)  +  l}  .  (22) 


For  the  quasi-harmonic  idealization,  defined  by  statements  5a  and  5b,  eqn  (22)  reduces  to 


E  V 

(pv^^y  =  -  ” 


2/3 


(M,)  {ln(^/^o)  ^  l}  • 


(23) 


If  the  tfr  functioais  quantified  as  ilr  =  (Fj^Vf^)  V■^^  in  accordance  with  the  definition  of  \|f  and 
statement  5b,  then  the  dF/dX  function  is  given  as 


dFIdk  =  3y2/3(p^y2/3y  ^ 


34f': 


(0/^^)  {ln(^/^(,)  +  l}  • 


(24) 


We  see  here  why  the  quasi-harmonic  case  is  an  important  one — ^namely,  because  the  interatomic 
force  and  its  spacial  derivatives  vary  only  with  the  characteristic  frequency  and  not  explicitly  as 
a  function  of  lattice  spacing  (X,  or  V).  The  relative  ratios  of  the  zero-temperature  vibrational  and 
volumetric  stiffnesses  may  be  compared  as  a  function  of  characteristic  temperature  ratio  and  are 
displayed  in  Figure  1 .  The  upper  curve,  for  vibrations,  is  parabolic  since  (according  to  statements 
lb  and  2)  it  is  the  spring  constant  which  governs  the  characteristic  frequency  (and  thus 
characteristic  temperature).  The  lower  curve,  derived  from  Segletes’  model  for  the  quasi¬ 
harmonic  case,  relates  the  relative  volumetric  stiffness  to  the  characteristic  temperature.  Were 
these  two  curves  to  overlay,  the  model  of  Dugdale  and  MacDonald  [8]  would  prevail.  As  it  is, 
the  vibrational  and  volumetric  stiffnesses  for  a  quasi-harmonic,  3-D  lattice  do  not  remain 
proportional,  though  they  do  so  approximately  near  the  reference  density. 

Eqn  (24)  expresses  the  volumetric  stiffness  in  terms  of  the  quasi-harmonic,  characteristic 
temperature  of  the  lattice.  A  converse  relationship  may  be  obtained  by  starting  with  eqn  (18), 
and  making  the  substitution  in  the  bracketed  term  that,  for  the  quasi-harmonic  case. 
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=  2/3  (tjr/V).  Employing  the  fact  that  F=p,V^\  the  resulting  equation  may  be  rearranged 
to  give 

-1  =  ^  (25) 

(tj/F'  +  F) 

The  left  side  of  the  equation  is,  by  way  of  eqn  (13),  (^7^ .  Similarly,  the  numerator  of  the  right 
side  of  the  equation  is  also  the  volume  derivative  of  the  denominator.  Thus,  we  see  that 

^  ~  tj/F'  +  F  =  —  +  F  ,  (26) 

3  f  ^ 

which  expresses  the  characteristic  frequency  in  terms  of  the  lattice  cell  force  and  its  spacial 
derivative. 

To  put  the  comparison  in  yet  another  perspective,  in  terms  of  the  moduli,  we  compute,  from 
eqns  (20)  and  (24),  that  Likewise,  from  eqn  (21)  and 

statements  lb  and  2,  we  know  that  J^/J^p  =  We  will  assume  now  the  initial 

condition,  that  S^p  =  B^p.  This  assumption  corresponds,  for  a  3-D  lattice,  to  the  situation  where 

(dFldk)y  equals  3{dFJdx\^.  If  this  assumption  should  subsequently  prove  wrong,  the  sole 

effect  will  be  that  eqns  (27)  and  (30),  to  follow,  will  be  off  by  a  constant  multiplier.  The  force- 
based  volumetric  and  vibrational  moduli,  for  the  quasi-harmonic  case,  may  thus  be  directly 
related  as 


= 


[ln(0/^o)  +  1] 


,) 


(27) 
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A  graph  of  this  function  is  shown  in  Figure  2,  as  a  function  of  the  quasi-harmonic,  relative 
characteristic  temperature,  This  figure  indicates  the  functional  relationship  between  the 

force-based  volumetric  and  vibrational  moduli  for  the  quasi-harmonic  case. 

Alternately,  this  relationship  may  be  expressed  as  a  function  of  specific  volume  ratio,  by 
converting  characteristic  temperature  to  specific  volume.  By  starting  with  the  quasi-harmonic 
function,  given  by 


f  =  (Of'o)  y"""  ’  (28) 

one  may  integrate  according  to  eqn  (13)  to  ascertain  the  logarithm  of  the  quasi-harmonic 
frequency  function: 


ln(^/^p)  =  3  fp  [1  -  (V/Vp)!'^]  .  (29) 

Interestingly,  the  quasi-harmonic  eqn  (29)  is  proportional  to  the  negative  of  the  a*  parameter 
employed  by  Rose  et  al.  [3],  as  the  nondimensional  lattice  parameter  in  their  universal  cold 
curve.  The  model  of  Rose  et  al.,  which  was  the  inspiration  for  Segletes’  equation  of 
state  [13,  14],  expresses  cold  energy  in  terms  of  compression  and  is  seen  from  eqn  (29)  to  relate 
closely  to  the  quasi-harmonic  case  idealized  from  Segletes’  model.  However,  unlike  Segletes’ 
equation  of  state,  Rose’s  model  does  not  include  thermal  pressure  effects  required  for  a  complete 
equation  of  state  and  has  no  provisions  for  Segletes’  general  case,  which  is  not  quasi-harmonic. 

Exponentiation  of  eqn  (29)  provides  the  characteristic  temperature  ratio.  The  quasi-harmonic 
characteristic  temperature  ratio  tells  exactly  why  the  real  world  is  not  quasi-harmonic — ^because, 
in  the  quasi-harmonic  case,  there  is  a  maximum  limiting  frequency  ratio  at  vanishingly  small 

volumes,  given  by  =  exp(3f'p).  Such  a  frequency  limit  would  translate  into  a  cold- 

compression  energy  limit,  as  well.  According  to  various  high-compression  theories  that  include 
electronic  effects,  finite  cold-compression  energies  under  infinite  compression  do  not  occur  in 
actual  materials.  Such  limiting  compressive  behavior,  though  absent  in  actual  materials. 
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characterizes  both  the  model  of  Rose  et  al.  [3],  as  well  as  the  quasi-harmonic  idealization.  A 
subsequent  report  will  address  how  actual  materials  may  be  described  by  Segletes’  theory,  by 
departing  from  the  tenets  of  the  quasi-harmonic  idealization.  In  the  report  to  appear,  the  author 
intends  to  show  that  the  data  support  the  notion  that,  as  lattice  spacing  increases,  the  behavior 
of  all  lattices  approaches  the  quasi-harmonic  idealization,  and  that  a  departure  from  this  ideal 
occurs  as  the  lattice  spacing  becomes  small  enough  to  permit  interatomic  repulsive  interactions. 

Substitution  may  be  made  into  eqn  (27),  in  order  to  express  the  volumetric  to  vibrational 
modulus  ratio,  in  terms  of  specific  volume,  as 


^  _  <1  +  3t,[l  -  iV/Voy^]} 

"  "  ‘  exp{3fjl  -  (y/V,r^]} 


(30) 


This  relationship  is  depicted  in  Figure  3  with  the  reference,  quasi-harmonic  Griineisen  parameter 
value,  fp,  as  a  parameter.  The  figure  shows  that,  only  as  fg  is  made  to  approach  zero,  does  the 

quasi-harmonic,  zero-temperature,  volumetric  to  vibrational  modulus  ratio  approach  a  constant 
value  (of  unity).  The  quasi-harmonic  lattice  force,  F=(p^V^^),  may  also  be  expressed  in  terms 
of  lattice  spacing,  rather  than  frequency,  through  the  use  of  eqns  (15)  and  (29)  and  the  definition 
that?i=F''': 


It  is  also  revealing  to  study  the  quasi-harmonic,  cold-energy  curve,  not  as  a  function  of 
frequency,  as  given  in  eqn  (16),  but  likewise  in  terms  of  relative  lattice  spacing,  using  eqn  (29) 
and  the  relation  between  X  and  V.  Figure  4  plots  this  relation  with  the  reference  quasi-harmonic 

Griineisen  parameter,  as  a  parameter.  The  energy  has  been  normalized  by  its  limiting  (V  =  0) 

cold-compression  energy,  which  is  finite  for  the  quasi-harmonic  case,  but  will  not  be  so  for  real 
materials.  The  figure  clearly  shows  that,  as  the  reference  Griineisen  parameter  approaches  zero, 
the  quasi-harmonic  idealization  of  Segletes’  equation  approaches  the  harmonic  case,  consistent 
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with  Dugdale  and  MacDonald  [8]  theory,  as  already  shown  for  T = 0.  Kittel  [23]  indicates  that, 
according  to  harmonic  theory,  there  is  no  thermal  expansion,  and  the  elastic  constants  are 
independent  of  pressure  (and  temperature).  A  zero  Griineisen  function  follows  from  a  zero 
thermal-expansion  coefficient.  Thus,  the  constant  ratio  of  the  elastic  moduli  in  Figure  3,  for  the 

limiting  case  of  =  0,  is  further  evidence  of  compatibility  with  harmonic  theory.  Furthermore, 

dF^ii,ldx  and  dF/dk  are  proportional  for  only  this  limiting  case  and  this  proportionality  is  exactly 
what  would  occur  for  a  3-D  lattice  of  atomic  masses  connected  by  simple  linear  springs.  In 

short,  the  quasi-harmonic  idealization  approaches  the  harmonic  limit  as  f  approaches  zero. 

The  harmonic  condition  has  been  traditionally  expressed  by  stipulating  the  potential  energy 
well  to  be  parabolic, 

d^E 

-Zfl  =  0  ,  (32) 

dk^ 


thus  resulting  in  a  constant  characteristic  frequency  and  an  identically  zero  Griineisen  function. 
We  see  now  that  the  harmonic  approximation  is  merely  the  limiting  case  of  the  quasi-harmonic 
idealization  of  Segletes’  equation  of  state,  when  the  value  of  the  Griineisen  function  is  made  to 
approach  zero  in  the  limit.  The  quasi-harmonic  potential,  given  by  eqn  (16),  satisfies  the 
following  relation: 


/72  E 

Cl  pot 


=  0  , 


(33) 


where  Epgp  the  potential  energy  weU  of  the  lattice,  is  the  lattice  cold  energy  relative  to  the  mfinite 
lattice  separation  condition,  or  Ep„  ==Ef.- E^. 
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5.  Conclusions 


The  historical  use  of  the  harmonic  approximation  in  the  study  of  lattice  dynamics  has  been 
briefly  reviewed  to  include  the  development  of  the  Dugdale-MacDonald  model  of  anharmonic 
lattice  behavior.  It  is  shown  how,  in  the  case  of  the  model  of  Dugdale  and  MacDonald  [8], 
results  rely  on  a  certain  assumption  about  the  governing  elasticity — ^an  assumption  that  follows 
from  1-D  theory  wherein  the  vibrational  and  volumetric  stiffnesses  are  proportional.  This 
assumption  has  been  widely  adopted  by  prior  treatments  of  the  subject,  but  has  been  shown  to 
be  inadequate  to  describe  the  behavior  of  a  3-D  lattice  at  large  volumetric  strain,  perhaps  as  a 
result  of  lateral  (Poisson  ratio)  coupling,  transverse  (shear)  waves  and/or  the  presence  of  non- 
nearest  neighbor  interatomic  interactions. 

The  relationship  between  vibrational  and  volumetric  moduli  has  been  studied  using  the  quasi- 
harmonic  idealization  of  Segletes’  equation  of  state.  The  quasi-harmonic  idealization  is  one  in 
which  the  form  of  Segletes’  equation  of  state  is  greatly  simplified,  while  the  Griineisen  function 
takes  on  a  specified,  idealized,  functional  behavior.  In  the  quasi-harmonic  idealization,  the 
interatomic  force  and  its  spacial  derivatives  vary  only  with  the  characteristic  frequency  of  the 
lattice  and  not  explicitly  with  the  lattice  spacing. 

Results  indicate  that  the  relationship  between  the  vibrational  and  volumetric  moduli  (and  their 
associated  natural  frequencies)  are  of  different  form.  In  the  vibrational  case,  the  stiffness, 
according  to  vibration  theory,  goes  as  (O^,  whereas  for  volumetric  distortions  of  a  3-D  lattice,  it 
may  be  inferred  from  the  idealized  form  of  Segletes’  equation  of  state  that  the  volumetric 
stiffness  goes  as  0)[1  +ln((o/(Oo)]-  It  is  this  disparity  that  causes  the  primary  difference  between 
the  model  of  Segletes  and  previous  work  by  Dugdale  and  MacDonald  [8]. 

The  idealized  Griineisen  relationship  that  governs  quasi-harmonic  behavior  is  such  that  a 
quasi-harmonic  lattice  reaches  a  limiting  cold-compression  energy  (and  associated  vibrational 
frequency)  in  the  limit  of  vanishingly  small  volume.  Such  idealized  behavior,  though  not 
indicative  of  real  materials,  is  compatible  with  the  theory  of  harmonic  lattice  vibrations.  This 
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report  explicitly  shows  that  the  quasi-harmonic  idealization,  for  the  limiting  case  where  the 
Griineisen  function  approaches  zero,  actually  approaches  the  harmonic  approximation. 

In  a  subsequent  report,  the  author  intends  to  address  how  real  material  behavior  may  be 
described  in  the  framework  of  Segletes’  equation  of  state,  by  departing  from  the  quasi-harmonic 
idealization,  and  how,  as  lattice  spacing  increases,  the  quasi-harmonic  idealization  is  approached 
for  real  lattices.  Such  a  result  is  wholely  compatible  with  the  conclusions  of  Plendl  [11],  who 
asserts  that  anharmonicity  is  “mainly  caused  by  the  forces  of  repulsion.,.,  but  is  close  to  zero 
when  ions  deviate  from  one  another,”  (i.e.,  when  X>Xo). 
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10.0 


Figure  1.  Relative  values  of  the  vibrational  (dF^,ji,ldx)  and  volumetric  (dFIdX)  stiffnesses,  as  a 
function  of  the  relative  characteristic  temperature  (or  frequency),  for  a  material 
obeying  the  quasi-harmonic  idealization. 


Figure  2.  Ratio  of  volumetric  to  vibrational  (force-based)  moduli,  as  a  function  of  the  relative 
characteristic  temperature  (or  frequency),  for  a  material  obeying  the  quasi-harmonic 
idealization. 
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Figure  3.  Ratio  of  volumetric  to  vibrational  (force-based)  moduli,  as  a  function  of  relative 
volume,  for  a  material  obeying  the  quasi-harmonic  idealization.  Six  curves  are 
shown,  corresponding  to  different  values  of  the  reference  Griineisen  parameter,  Fq. 


X/Xo 


Figure  4.  Relative  cold  compression  energy  versus  relative  lattice  spacing,  with  reference 
Griineisen  value,  Fq,  as  a  parameter,  for  quasi-harmonic  material.  Fq  approaches 
zero,  harmonicity  is  approached. 
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ATTN  KEZHUN  LI 
BERKELEY  CA  94720 

1  UC  DAVIS 

INST  OF  THEORETICAL  DYNAMICS 
ATTN  E  G  PUCKETT 
DAVIS  CA  95616 

1  UC  LOS  ANGELES 

DEPT  OF  MAT  SCIENCE  &  ENGNG 
ATTN  J  J  GILMAN 
LOS  ANGELES  CA  90024 
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ORGANIZATION 

COPIES 
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2 

UC  SAN  DIEGO 

DEPT  APPL  NECH  &  ENGR 

SVCS  ROll 

ATTN  S  NEMAT-NASSER 

M  MEYERS 

1 

UNIVERSITY  OF  TEXAS 

DEPT  OF  MECHANICAL  ENGINEERING 
ATTN  ERIC  P  FAHRENTHOLD 

AUSTIN  TX  78712 

LA  JOLLA  CA  92093-041 1 

1 

VIRGINIA  POLYTECHNIC  INSTITUTE 
COLLEGE  OF  ENGINEERING 

2 

UNIV  OF  ALA  HUNTSVILLE 

AEROPHYSICS  RSCH  CTR 

ATTN  GARY  HOUGH 

ATTN  R  BATRA 

BLACKSBURG  VA  24061-0219 

DAVID  J  LIQUORNIK 

PO  BOX  999 

HUNTSVE:.LE  AL  35899 

2 

AEROJET 

ATTN  J  CARLEONE 

S  KEY 

PO  BOX  13222 

1 

UNIV  OF  ALA  HUNTSVILLE 

CIVIL  ENGRNG  DEPT 

SACRAMENTO  CA  95813-6000 

ATTN  WILLIAM  P  SCHONBERG 
HUNTSVILLE  AL  35899 

2 

AEROJET  ORDNANCE 

ATTN  P  WOLF 

G  PADGETT 

1 

UNIVERSITY  OF  CHICAGO 

DEPT  OF  THE  GEOPHYSICAL  SCIENCES 
ATTN  G  H  MILLER 

1100  BULLOCH  BLVD 

SOCORRO  NM  87801 

5734  S  ELLIS  AVE 

CHICAGO  IL  60637 

3 

ALLIANT  TECHSYSTEMS  INC 

ATTN  T  HOLMQUIST  MNl  1  2720 

R  STRYK 

3 

UNIVERSITY  OF  DAYTON  RSCH  INST 
KLA14 

ATTN  N  BRAR 

D GROVE 

G  R  JOHNSON  MNll  2925 

600  SECOND  ST  NE 

HOPKINS  MN  55343 

4 

A  PIEKUTOWSKI 

300  COLLEGE  PARK 

DAYTON  OH  45469-0182 

UNIVERSITY  OF  DELAWARE 

1 

ALME  AND  ASSOCIATES 

ATTN  MARVIN  L  ALME 

6219  BRIGHT  PLUME 

COLUMBIA  MD  21044-3790 

DEPT  OF  MECHANICAL  ENGINEERING 
ATTN  PROF  J  GILLESPIE 

DEAN  R  B  PIPES 

PROF  J  VINSON 

PROF  D  WILKINS 

1 

APPLIED  RESEARCH  ASSOC  INC 

ATTN  JEROME  D  YATTEAU 

5941  S  MIDDLEFIELD  RD  SUITE  100 
LITTLETON  CO  80123 

NEWARK  DE  19716 

2 

APPLIED  RESEARCH  ASSOC  INC 

ATTN  DENNIS  GRADY 

1 

UNIVERSITY  OF  MARYLAND 

ATTN  PHYSICS  DEPT  (BLDG  082) 

COLLEGE  PARK  MD  20742 

FRANK  MAESTAS 

4300  SAN  MATEO  BLVD  SE 
ALBUQUERQUE  NM  87110 

1 

UNIVERSITY  OF  PUERTO  RICO 

DEPT  CHEMICAL  ENGINEERING 

ATTN  L  A  ESTEVEZ 

MAYAGUEZ  PR  00681-5000 

1 

BATTELLE 

ATTN  ROBER  M  DUGAS 

7501  S  MEMORIAL  PKWY  SUITE  101 
HUNTSVILLE  AL  35802-2258 
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3  BOEING  AEROSPACE  CO 

SHOCK  PHYSICS  &  APPLIED  MATH 
ENGINEERING  TECHNOLOGY 
ATTN  R  HELZER 
T  MURRAY 
J  SHRADER 
PO  BOX  3999 
SEATTLE  WA  98124 

1  BOEING  HOUSTON  SPACE  STN 
ATTN  RUSSELL  F  GRAVES 
BOX  58747 
HOUSTON  TX  77258 

1  BRIGS  CO 

ATTN  JOSEPH  E  BACKOFEN 
2668  PETERSBOROUGH  ST 
HERNDON  VA  20171-2443 

1  CALIFORNIA  RSCH  &  TECHNOLOGY 
ATTN  M  MAJERUS 
PO  BOX  2229 
PRINCETON  NJ  08543  . 

1  CENTURY  DYNAMICS  INC 
ATTN  N  BIRNBAUM 
2333  SAN  RAMON  VALLEY  BLVD 
SAN  RAMON  CA  94583-1613 

1  COMPUTATIONAL  MECHANICS 
CONSULTANTS 
ATTN  J  A  ZUKAS 
POBOX  11314 
BALTIMORE  MD  21239-0314 

1  CYPRESS  INTERNATIONAL 
ATTN  A  CAPONECCHI 
1201  E  ABINGDON  DR 
ALEXANDRIA  VA  22314 

1  DEFENSE  TECHNOLOGY  INTL.  INC 
ATTN  D  E  AYER 
THE  STARK  HOUSE 
22  CONCORD  ST 
NASHUA  NH  03060 

1  DESKIN  RESEARCH  GROUP  INC 
ATTN  EDWARD  COLLINS 
2270  AGNEW  RD 
SANTA  CLARA  CA  95054 


1  DOW  CHEMICAL  INC 
ORDNANCE  SYSTEMS 
ATTN  C  HANEY 
A  HART 
B  RAFANIELLO 
800  BUILDING 
MIDLAND  MI  48667 

1  G  E  DUVALL 

5814  NE  82ND  COURT 
VANCOUVER  WA  98662-5944 

3  DYNA  EAST  CORP 
ATTN  P  C  CHOU 
R  CICCARELLI 
WFLIS 

3620  HORIZON  DRIVE 
KING  OF  PRUSSIA  PA  19406 

3  DYNASEN 

ATTN  JACQUES  CHAREST 
MICHAEL  CHAREST 
MARTIN  LILLY 
20  ARNOLD  PL 
GOLETA  CA  93117 

1  R  J  EICHELBERGER 
409  W  CATHERINE  ST 
BEL  AIR  MD  21014-3613 

1  ELORET  INSTITUTE 

ATTN  DAVID  W  BOGDANOFF  MS  230  2 
NASA  AMES  RESEARCH  CENTER 
MOFFETT  FIELD  CA  94035 

3  ENIG  ASSOCIATES  INC 
ATTN  J  ENIG 
D  J  PASTINE 
M  COWPERTHWAITE 
SUITE  500 

11120  NEW  HAMPSHIRE  AVE 
SILVER  SPRING  MD  20904-2633 

1  EXPLOSIVE  TECHNOLOGY 
ATTN  M  L  KNAEBEL 
PO  BOX  KK 
FAIRFIELD  CA  94533 
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1  GB  TECH  LOCKHEED 
ATTN  JAY  LAUGHMAN 
2200  SPACE  PARK  SUITE  400 
HOUSTON  TX  77258 

2  GB  TECH  LOCKHEED 

ATTN  LUCILLE  BORREGO  C23C 
JOE  FALCON  JR  C23C 
2400  NASA  ROAD  I 
HOUSTON  TX  77058 

6  GDLS 

38500  MOUND  RD 
ATTN  W  BURKE  MZ436-2I-24 
G  CAMPBELL  MZ436-30-44 
D  DEBUSSCHER  MZ436-20-29 
J  ERIDON  MZ436-2I-24 
W  HERMAN  MZ  435-01-24 
S  PENTESCU  MZ436-21-24 
STERLING  HTS  MI  48310-3200 

2  GENERAL  RESEARCH  CORP 
ATTN  A  CHARTERS 
TMENNA 
PO  BOX  6770 

SANTA  BARBARA  CA  93160-6770 

2  GRC  INTERNATIONAL 

ATTN  TIMOTHY  M  CUNNINGHAM 
WILLIAM  M  ISBELL 
5383  HOLLISTER  AVE 
SANTA  BARBARA  CA  931 1 1 

6  INST  OF  ADVANCED  TECHNOLOGY 
UNIVERSITY  OF  TX  AUSTIN 
ATTN  STEPHEN  J  BLESS 
JAMES  CAZAMIAS 
HARRY  D  FAIR 
THOMAS  M  KIEHNE 
DAVID  LITTLEFIELD 
MIKE  NORMANDIA 
4030-2  W  BRAKER  LN 
AUSTIN  TX  78759 

1  INTERNATIONAL  RESEARCH  ASSOC 
ATTN  D  ORPHAL 
4450  BLACK  AVE 
PLEASANTON  CA  94566 


NO.  OF 

COPIES  ORGANIZATION 

1  INTERPLAY 

ATTN  F  E  WALKER 
18  SHADOW  OAK  RD 
DANVILLE  CA  94526 

1  KAMAN  SCIENCES  CORP 
ATTN  DENNIS  L  JONES 
2560  HUNTINGTON  AVE  SUITE  200 
ALEXANDRIA  VA  22303 

8  KAMAN  SCIENCES  CORP 
ATTN  J  ELDER 
RICHARD  P  HENDERSON 
DAVID  A  PYLES 
FRANK  R  SAVAGE 
JAMES  A  SUMMERS 
JAMES  S  WILBECK 
TIMOTHY  W  MOORE 
THY  YEM 

600  BLVD  S  SUITE  208 
HUNTSVILLE  AL  35802 

3  KAMAN  SCIENCES  CORP 
ATTN  SHELDON  JONES 
GARY  L  PADEREWSKI 
ROBERT  G  PONZINI 

1500  GRDN  OF  THE  GODS  RD 
COLORADO  SPRINGS  CO  80907 

4  KAMAN  SCIENCES  CORP 
ATTN  NASrr  ARI 
STEVE  R  DIEHL 
WILLIAM  DOANE 
VERNON  M  SMITH 

PO  BOX  7463 

COLORADO  SPRINGS  CO  80933-7463 

1  D  R  KENNEDY  &  ASSOC  INC 
ATTN  D  KENNEDY 
PO  BOX  4003 

MOUNTAIN  VIEW  CA  94040 

1  KERLEY  PUBLISHING  SERVICES 
ATTN  G  I  KERLEY 
PO  BOX  13835 

ALBUQUERQUE  NM  87192-3835 


35 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


2  KTECH  CORPORATION 
ATTN  FRANK  W  DAVIES 
LARRY  M  LEE 

901  PENNSYLVANIA  NE 
ALBUQUERQUE  NM  871 10 

I  LIVERMORE  SOFTWARE  TECH  CORP 
ATTN  J  O  HALLQUIST 
2876  WAVERLY  WAY 
LIVERMORE  CA  94550 

I  LOCKHEED  MARTIN  MISSLE  &  SPACE 
ATTN  WILLIAM  R  EBERLE 
PO  BOX  070017 
HUNTSVILLE  AL  35807 

3  LOCKHEED  MARTIN  MISSILE  &  SPACE 
ATTN  M  A  LEVIN  ORG  81  06  BLDG  598 
M  R  MCHENRY 

T  A  NGO  ORG  81  10  BLDG  157 
111  LOCKHEED  WAY 
SUNNYVALE  CA  94088 

4  LOCKHEED  MISSILE  &  SPACE  CO 
ATTN  JOHN  R  ANDERSON 
WILLIAM  C  KNUDSON 

S  KUSUMI  0  81  11  BLDG  157 
J  PHILLIPS  0  54  50 
PO  BOX  3504 
SUNNYVALE  CA  94088 

1  LOCKHEED  MISSILE  &  SPACE  CO 
ATTN  R  HOFFMAN 
SANTA  CRUZ  FACILITY 
EMPIRE  GRADE  RD 
SANTA  CRUZ  CA  95060 

1  LOCKHEED  NASA  JSC 
SPACE  SCIENCE  BRANCH 
ATTN  JAMES  HYDE 
BOX  58561  MC  B22 
HOUSTON  TX  77258 

1  LOCKHEED  MARTIN  AEROSPACE 
ATTN  D  R  BRAGG 
PO  BOX  5837  MP  109 
ORLANDO  FL  32855 


1  MCDONNELL  DOUGLAS 
ASTRONAUTICS  CO 
ATTN  B  L  COOPER 
5301  BOLSA  AVE 
HUNTINGTON  BEACH  CA  92647 

1  ORLANDO  TECHNOLOGY  INC 
ATTN  DANIEL  A  MATUSKA 
PO  BOX  855 
SHALIMAR  FL  32579 

1  PHYSICAL  SCIENCES  INC 
ATTN  PETER  NEBOLSINE 
20  NEW  ENGLAND  BUS  CTR 
ANDOWER  MA  01810 

3  PHYSICS  INTERNATIONAL 
ATTN  R  FUNSTON 
G  FRAZIER 
L  GARNETT 
PO  BOX  5010 
SAN  LEANDRO  CA  94577 

1  PRC  INC 

ATTN  J  ADAMS 

5166  POTOMAC  DR  #103 

KING  GEORGE  VA  22485-5824 

1  RAYTHEON  ELECTRONIC  SYSTEMS 
ATTNR  KARPP 
50  APPLE  HELL  DRIVE 
TEWKSBURY  MA  01876 

1  ROCKWELL  INTERNATIONAL 
ROCKETDYNE  DIVISION 
ATTN  H  LEIFER 
16557  PARK  LN  CIRCLE 
LOS  ANGELES  CA  90049 

1  ROCKWELL  MISSILE  SYS  DIV 
ATTN  T  NEUHART 

1800  SATELLITE  BLVD 
DULUTH  GA  30136 

2  SAIC 

ATTN  JAMES  FURLONG 
GREGORY  J  STRAUCH 
1710  GOODRIDGE  DR 
MCLEAN  VA  22102 
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1  SAIC 

ATTN  MICHAEL  W  MCKAY 
10260  CAMPUS  POINT  DR 
SAN  DIEGO  CA  92121 

1  SHOCK  TRANSIENTS  INC 
ATTN  DAVID  DAVISON 
BOX  5357 

HOPKINS  MN  55343 

2  SIMULATION  &  ENG  CO  INC 
ATTN  ELSA  I  MULLINS 
STEVEN  E  MULLINS 

8840  HWY  20  SUITE  200  N 
MADISON  AL  35758 

2  SOUTHERN  RESEARCH  INSTITUTE 
ATTN  LINDSEY  A  DECKARD 
DONALD  P  SEGERS 
PO  BOX  55305 

BIRMINGHAM  AL  35255-5305 

5  SRI  INTERNATIONAL 
ATTN  JAMES  D  COLTON 
D CURRAN 
R  KLOOP 
R  L  SEAMAN 
D  A  SHOCKEY 
333  RAVENSWOOD  AVE 
MENLO  PARK  CA  94025 

2  TELEDYNE  BROWN  ENGR 
ATTN  JIM  W  BOOTH 
MARTIN  B  RICHARDSON 
PO  BOX  070007  MS  50 
HUNTSVILLE  AL  35807-7007 

1  ZERNOW  TECHNICAL  SVCS  INC 
ATTN  LOUIS  ZERNOW 
425  W  BONITA  AVE  SUITE  208 
SAN  DIMAS  CA  91773 


NO.  OF 

COPIES  ORGANIZATION 
ABERDEEN  PROVING  GROUND 

30  DIR,  USARL 

ATTN:  AMSRL-WM,  I  MAY 
AMSRL-WM-MF, 

S  CHOU 
D  DANDEKAR 
A  RAJENDRAN 
AMSRL-WM-PB,  A  ZIELINSKI 
AMSRL-WM-PD,  G  GAZONAS 
AMSRL-WM-T,  W  F  MORRISON 
AMSRL-WM-TA, 

M  BURKINS 
W  GILLICH 
W  BRUCHEY 
J  DEHN 
G  HLBEY 
W  A  GOOCH 
H  W  MEYER 
E  J  RAPACKI 
J  RUNYEON 
AMSRL-WM-TB, 

R  FREY 
P  BAKER 
J  STARKENBERG 
AMSRL-WM-TC, 

W  S  DE  ROSSET 
T  W  BJERKE 
R  COATES 
F  GRACE 
K  KIMSEY 
M  LAMPSON 
D  SCHEFFLER 
S  SCHRAML 
B  SORENSEN 
R  SUMMERS 
W  WALTERS 
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15  DIR,  USARL 

AMSRL-WM-TD, 

A  M  DIETRICH 
K  FRANK 
J  HARRISON 
M  RAFTENBERG 
G  RANDERS-PEHRSON 
M  SCHEIDLER 
S  SCHOENFELD 
S  SEGLETES  (5  CP) 

J  WALTER 
T  WRIGHT 

AMSRL-WM-WD,  A  PRAKASH 


NO.  OF 
COPIES 

3 


1 


1 


1 


1 


1 


1 


ORGANIZATION 

AERONAUTICAL  &  MARITIME 
RESEARCH  LABORATORY 
ATTN  N  BURMAN 
S  CIMPOERU 
D  PAUL 
PO  BOX  4331 
MELBOURNE  VIC  3001 
AUSTRALIA 

EMBASSY  OF  AUSTRALIA 
ATTN  R  WOODWARD 
COUNSELLOR  DEFENCE  SCIENCE 
1601  MASSACHUSETTS  AVE  NW 
WASHINGTON  DC  20036-2273 

ABTEILUNG  FUER  PHYSIKALISCHE 
CHEMIE 

MONTANUNIVERSITAET 
ATTN  E  KOENIGSBERGER 
A  8700  LEOBEN 
AUSTRIA 

PRB  S  A 

ATTN  M  VANSNICK 

AVENUE  DE  TERVUEREN  168  BTE  7 

BRUSSELS  B  1150 

BELGIUM 

ROYAL  MILITARY  ACADEMY 
ATTN  E  CELENS 
RENAISSANCE  AVE  30 
B1040  BRUSSELS 
BELGIUM 

BULGARIAN  ACADEMY  OF  SCIENCES 
SPACE  RESEARCH  INSTITUTE 
ATTN  VALENTIN  GOSPODINOV 
1000  SOFIA  PO  BOX  799 
BULGARIA 

CANADIAN  ARSENALS  LTD 
ATTN  P  PELLETIER 
5  MONTEE  DES  ARSENAUX 
vn  .T  .IF.  DE  GRADEUR  PQ  J5Z2 
CANADA 

DEFENCE  RSCH  ESTAB  SUFFIELD 
ATTN  D  MACKAY 

RALSTON  ALBERTA  TOJ  2NO  RALSTON 
CANADA 


NO.  OF 

COPIES  ORGANIZATION 

1  DEFENCE  RSCH  ESTAB  SUFFIELD 
ATTN  CHRIS  WEICKERT 
BOX  4000  MEDICINE  HAT 
ALBERTA  TIA  8K6 
CANADA 

1  DEFENCE  RSCH  ESTAB  VALCARTIER 
ARMAMENTS  DIVISION 
ATTN  R  DELAGRAVE 
2459  PIE  XI  BLVD  N 
PO  BOX  8800 

CORCELETTE  QUEBEC  GOA  IRO 
CANADA 

1  UNIVERSITY  OF  GUELPH 
PHYSICS  DEPT 
ATTN  C  G  GRAY 
GUELPH  ONTARIO 
NIG  2W1 
CANADA 

1  CEA 

ATTN  ROGER  CHERET 
CEDEX  15 

313  33  RUE  DE  LA  FEDERATION 

PARIS  75752 

FRANCE 

1  CEA 

CISI  BRANCH 
ATTN  PATRICK  DAVID 
CENTRE  DE  SACLAY  BP  28 
GIF  SUR  YVETTE  91192 
FRANCE 

1  CEA/CESTA 

ATTN  ALAIN  GEILLE 
BOX  2  LE  BARP  33114 
FRANCE 

6  CENTRE  D’ETUDES  DE  GRAMAT 
ATTN  SOLVE  GERARD 
CHRISTIAN  LOUPIAS 
PASCALE  OUTREBON 
J  CAGNOUX 
C  GALLIC 
J  TRANCHET 
GRAMAT  46500 
FRANCE 
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2 


3 


6 


1 


1 


ORGANIZATION 
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CENTRE  D’ETUDES  DE  LIMEIL-VALENTON 
ATTN  CHRISTIAN  AUSSOURD 
JEAN-CLAUDE  BOZIER 
SAINT  GEORGES  CEDEX 
VILLENEUVE  94195 
FRANCE 

CENTRE  D’ETUDES  DE  VAUJOURS 

ATTN  PLOTARD  JEAN-PAUL 

ERIC  BOTTET 

TAT  SIHN  VONG 

BOITE  POSTALE  NO  7 

COUNTRY  77181 

FRANCE 

CENTRE  DE  RECHERCHES 
ET  D’ETUDES  D’ARCUEIL 
ATTN  D  BOUVART 
C  COTTENNOT 
S  JONNEAUX 
HORSINI 
S  SERROR 
FTARDIVAL 

16  BIS  AVENUE  PRIEUR  DE 
LA  COTE  D’OR 
F94114  ARCUEIL  CEDEX 
FRANCE 

DAT  ETBS  CETAM 

ATTN  CLAUDE  ALTMAYER 

ROUTE  DE  GUERRY  BOURGES 

18015 

FRANCE 

ETBS  DSTI 
ATTN  P  BARNIER 
ROUTE  DE  GUERAY 
BOITE  POSTALE  712 
18015  BOURGES  CEDEX 
FRANCE 

FRENCH  GERMAN  RESEARCH  INST 
ATTN  CHANTERET  P-Y 
CEDEX  12  RUE  DE  I’lNDUSTRIE 
BP  301 

F68301  SAINT-LOUIS 
FRANCE 


5  FRENCH  GERMAN  RESEARCH  INST 
ATTN  HANS-JURGEN  ERNST 
FRANCIS  JAMET 
PASCALE  LEHMANN 
K  HOOG 
H  LERR 

CEDEX  5  5  RUE  DU  GENERAL 
CASSAGNOU 
SAINT  LOUIS  68301 
FRANCE 

1  LABORATOIRE  DE  TECHNOLOGIE  DES 
SURFACES 

ECOLE  CENTRALE  DE  LYON 
ATTN  VINET  P 
BP  163 

69131  ECULLY  CEDEX 
FRANCE 

1  BATTELLE  INGENIEUTECHNIK  GMBH 
ATTN  W  FUCHE 
DUESSELDORFFER  STR  9 
ESCHBORN  D  65760 
GERMANY 

1  CONDAT 

ATTN  J  KIERMEIR 
MAXIMILIANSTR  28 
8069  SCHEYERN  FERNHAG 
GERMANY 

1  DEUTSCHE  AEROSPACE  AG 
ATTN  MANFRED  HELD 
POSTFACH  13  40 
D  86523  SCHROBENHAUSEN 
GERMANY 

1  DIEHL  GBMH  AND  CO 
ATTN  M  SCHELDKNECHT 
nSCHBACHSTRASSE  16 
D  90552  ROTBENBACH  AD  PEGNITZ 
GERMANY 
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COPIES  ORGANIZATION 

5  ERNST  MACH  INSTITUT 
ATTN  VOLKER  HOHLER 
E  SCHMOLINSKE 
E  SCHNEIDER 
A  STILP 
KTHOMA 
ECKERSTRASSE  4 
D-7800  FREIBURG  I  BR  791  4 
GERMANY 

1  EUROPEAN  SPACE  AGENCY 
ATTN  WALTER  FLURY 
ESOC  5 

ROBT  BOSCHT  STRASSE 
DARMSTADT  6100 
GERMANY 

3  FRAUNHOFER  INSTITUT  FUER 
KURZZEITDYNAMIK 
ERNST  MACH  INSTITUT 
ATTN  H  ROTHENHAEUSLER 
H  SENF 

E  STRASSBURGER 
HAUPTSTRASSE  18 
D79576  WEIL  AM  RHEIN 
GERMANY 

3  FRENCH  GERMAN  RESEARCH  INST 
ATTN  HARTMUTH  F  LEHR 
ROLF  HUNKLER 
ERICH  WOLLMANN 
POSTFACH  1260 
WEIL  AM  RHEIN  D-79574 
GERMANY 

2  lABG 

ATTN  M  BORRMANN 
H  G  DORSCH 
EINSTEINSTRASSE  20 
D  8012  OTTOBRUN  B  MUENCHEN 
GERMANY 

1  INGENIEURBURO  DEISENROTH 
AUF  DE  HARDT  33  35 
D5204  LOHMAR  1 
GERMANY 


NO.  OF 

COPIES  ORGANIZATION 

3  TU  CHEMNITZ-ZWICKAU 
ATTN  I  FABER 
L  KRUEGER 
LOTHAR  MEYER 

FAKULTAET  FUER  MASCHINENBAU  U. 

VERFAHRENSTECHNIK 
SCHEFFELSTRASSE  110 
09120  CHEMNITZ 
GERMANY 

1  TECHNISCHE  UNIVERSITAT  MUENCHEN 
ATTN  EDUARD  B  IGENBERGS 
RICHARD  WAGNER  STR  18  111 
MUENCHEN  2  D8000 
GERMANY 

1  BHABHA  ATOMIC  RESEARCH  CENTRE 
HIGH  PRESSURE  PHYSICS  DIVISION 
ATTN  N  SURESH 
TROMBAY  BOMBAY  400  085 
INDIA 

5  RAFAEL  BALLISTICS  CENTER 
ATTN  EREZ  DEKEL 
YEHUDA  PARTOM 
G  ROSENBERG 
Z  ROSENBERG 
Y  YESHURUN 
PO  BOX  2250 
HAIFA  31021 
ISRAEL 

1  TECHNION  INST  OF  TECH 
FACULTY  OF  MECH  ENGNG 
ATTN  SOL  BODNER 
TECHNION  CITY 
HAIFA  32000 
ISRAEL 

1  CRC  RESEARCH  INST  INC 

ATTN  MASAHIDE  KATAYAMA 
STRUCTURAL  ENGINEERING  DEPT 
1  3  316  NAKASE  MIHAMA  KU 
CHIBA  SHI  261  01 
JAPAN 
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